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ABSTRACT 

We present a detailed study of the surface brightness profiles of dense filaments in IC 5146 using recent Herschel 
observations done with SPIRE. We describe the profile through an equilibrium solution of a self-gravitating isothermal 
cylinder pressure confined by its surrounding medium. In this first analysis we applied a simple modified black body 
function for the emissivity, neglecting any radiative transfer effects. Overall we found a good agreement of the observed 
surface brightness profiles with the model. The filaments indicate strong self-gravity with mass line densities M/l >~ 
0.5(M/Z) max where (M/Z) max is the maximum possible mass line density. In accordance with the model expectations we 
found a systematic decrease of the FWHM, a steepening of the density profile, and for filaments heated by the interstellar 
radiation field a decrease of the luminosity to mass ratio for higher central column density and mass line density. We 
illustrate and discuss the possibility of estimating the distance, external pressure, and dust opacity. For a cloud distance 
D ~ 500 pc and a gas temperature of T cy \ — 10 K the model implies an external pressure p cx t/fc ~ 2 x 10 4 K cm" 3 and 
an effective dust emission coefficient at 250 /im given by Sk™ ~ 0.0588 cm 2 g" 1 where 5 is the dust-to-gas ratio. Given 
the largest estimate of the distance to the cloud complex, 1 kpc, the model yields an upper limit Sk™ ~ 0.12 cm 2 g _1 . 

Key words. Methods: observational, Techniques: photometric, ISM: clouds, ISM: dust, extinction, ISM: structure, ISM: 
individual objects: IC 5146, Submillimeter: ISM 



1. Introduction 

Dense self-gravitating structures are an essential part of 
the interstellar medium and are the very first steps in the 
star formation process. Cold self-gravitating structures in 
the interstellar medium (ISM) can be found in spherical, 
elliptical or elongated for ms. Cold spherical fo rms are re- 
ferred to as Bok Globules (|Bok fc Reillvlll947l> if they ap- 
pear rather isolated in t he ISM. They are, as pointed out by 
ICurrv fc McKed (|2000l ). physically similar to the so-called 
condensed cores found in molecular clouds. A still larger 
fraction of the cold gas seems to be located in filamentary 
structures which are found to be a common phenomena 
throughout the ISM. They are seen in non star forming 
clouds such as Polaris as well a s in star forming cloud such 
as Aquila (I Andre et all 120101: IMen'shchikov et al.1 l2010t 
lArzoumanian et all 120111 ). Like the Bok Globules or cores 
they are well defined dense structures which make them 
perfect objects to study the cold molecular phase. 

As most of the star forming cores are found in the fil- 
amentary structures it is discussed whether there exists a 
strong link between the two which is important for the star 
formation process and responsible for the initial mass func- 
tion of the stars (IMF). In a common scenario the con- 
densed structures form hierarchically out of gravitational 
instabiliti es: the filaments out of sheets and the cores out of 
filaments (ISchneider fc Elmegreenlll979| iGaida et al.lll984t 



[19871: llnutsuka fc Mivamal Il992l [l997t ICurrv fc McKed 



2000: iFieee fc Pudritd boOOat IKandori et all 120051: 



Fischera fc Dopital 120081: iFischeral 120111) In our last 



Hanawa et all 119931: ICurrv fe McKed 120001: lAndre et all 
2010t IMen'shchikov et al.ll2010D . 

Cold clouds located in the interstellar medium 
are p r essure - confined by t he surr o undin g medium 
(|Ebertl 119551 : H onno 3 119561: iMcCreal 119571: iNagasawal 



paper (jFischera fc Martini ( 2012D . baper fl ) we have an- 
alyzed in detail the physical properties of pressurized 
cylinders in direct comparison with pressurized spheres. 

Filamentary structures in the cloud complexes IC 5146, 
Polaris, and Aquila show a typical FWHM ~ 0.1 pc but 
a large variation of t he central column density over three 
orders of magnitude (|Arzoumanian et al.l l201lh . The fila- 
ments in non star forming regions like Polaris appear to be 
smaller than the ones in star forming complexes as Aquila. 
The filaments of IC 5146 cover a physical regime in between 
the two other samples. The data overall indicate a trend 
that filaments are larger in size at higher column density. 

In |paper~H we have argued that the overall trend is prob- 
ably related to massive elongated structures with embed- 
ded dense filaments which cannot be described through a 
simple isothermal cloud model. However, the physical pa- 
rameters of elongated structures below a certain column 
density can be well explained by a model of an isothermal 
self-gravitating cylinder which is pressurized by the sur- 
rounding medium. 

This paper presents a pilo t study in which we apply 
the physical model outlined in Ipaper J to analyze the den- 
sity structure. The filaments chosen for study are in the 
IC 5146 field (GLON = 94 °, GLAT = -5°) studied 
by lArzoumanian et"al (12011) as part o f the Gould Belt 
Survey key project ( Andre et al. 2010) of the Herschel 
Space Observatory. IC 5146 is a key Herschel field in re- 
vealing the possible 'common' properties of filaments in 
nearby molecular clouds and therefore it is an important 
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one for which to develop physical models. We demonstrate 
how sub-millimeter observations can be used to derive an 
independent measurement of the cloud distance and of the 
ambient confining pressure. We further show that the model 
can be used to study the dust emission behavior in the sub- 
millimeter regime. _ ^ 

In a comparable study iKandori et al.l (|2005|) analyzed 
the column density profile of an ensemble of Bok-Globules. 
The authors found that a model of a self-gravitating pres- 
surized isothermal sphere can reproduce the observed col- 
umn density profiles which they derived using extinction 
measurements. They demonstrated that the model can be 
used to derive distances and pressures of the surrounding 
medium. 

In Sect. [5] we present the observational parameters for 
the filaments and the method used to extract surface bright- 
ness profiles. The physical model of the filaments is de- 
scribed in Sect. [3] The technique used to derive the physical 
parameters is described in Sect. 13.41 and the results for the 
individual filaments are presented in Sect.|U We discuss the 
results of the profile fitting and the estimates of the distance 
and pressure in Sect. [5] where we also address systematic 
uncertainties of the applied method. A short summary is 
given in Sect. [5] 

2. Observations 

The study is based on data from the Herschel Science 
Archived The Gould Belt Survey was mad e in parallel 
mode, with both P ACS dPoglitsch et al.l[2010h and SPIRE 
(jGriffin et al.l l2010h simultaneously, with two scans of the 
field in nearly orthogonal directions. Images made combin- 
ing the two scans self-consistently can lead to an optimal 
map (a 'level 2.5' product), minimizing the striping that 
appears, especially for PACS, in images made from a sin- 
gle scan direction (a 'level 2' product). As we began this 
study we were able to obtain the desired 'level 2.5' images 
for SPIRE, but even at the time of writing (late July 2012) 
these are not available for PACS data. We therefore decided 
to proceed with SPIRE data only. 

In the future the quality of archival PACS images 
will certainly improve (e.g., the promised 'Level 2.5' maps 
combing the cross scans, plus convergence results from dif- 
ferent map- making algorithmns) , which will in principle al- 
low additional constraints on the parmeters of the model. 
However, it is not necessarily simple. The dust emission 
depends more strongly on the temperature gradient in the 
filament at shorter wavelengths. Thus, although the reso- 
lution is higher at shorter wavelengths, the PACS data are 
more strongly affected by temperature variations which is 
actually expected to le ad to broader, not narrower, sur- 
face brightness profiles dFischera fc Dopital 120081 : iFischeral 
I201H iFischera fc Martini I2012D . This would make it nec- 
essary to use a more complex model. The profiles at the 
submm SPIRE bands are only slightly affected by temper- 
ature variations as the surface brightness in the Rayleigh- 
Jeans regime is proportional to dust temperature, not non- 
linearly to some higher power. This makes the SPIRE data 
used optimal for analyzing the structure using the simple 
single dust temperature model adopted here. 

The SPIRE instrument provides photometric images at 
reference wavelengths A r = 250, 350, and 500 [im. We 
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demonstrate that the high accuracy of the SPIRE data al- 
low a first analysis of the density structure using a sim- 
plified model and, despite the incomplete sampling of the 
spectral energy distribution (SED), even an estimate of the 
dust temperature and dust masses. 

The IC 5146 cloud complex has as the most prominent 
features a star formation region to the far east called the 
Cocoon Nebula and a massive elongated structure to the 
west known as Northern Streamer. The distance to the 
complex is not well establishe d with estimates varying by a 
factor of two from ~ 460 pc (|Lada et al.lll999j) to ~ 1 kpc 
(lHarvev et al.ll2008l). If we assume for the Sun a height 
zq = 26 ± 3 pc (jMaiaess et al.ll2009h above the mid-plane 
of the Galactic disc the position of the cloud complex lies 
somewhere between z — —14 pc and z = —61 pc below the 
mid-plane. 

lArzoumanian et al.l ([20111) combined SPIRE and PACS 
observations to produce column density maps for the same 
field which were then used to study the physical properties 
of individual filamentary structures including the column 
density profiles. In their analysis an automatic routine was 
applied to identify individual filamentary structures. The 
profiles were derived as an average of numerous orthogonal 
cuts through the filaments. 

We have chosen a different and a more traditional ap- 
proach to derive surface brightness profiles by estimating 
the profiles on the basis of a rectangular aperture laid over 
a section of a filament with well defined structure and back- 
ground. We have chosen four filaments named 'a', 'b', 'c', 
and 'd' in the paper. The filaments and the apertures are 
shown in Fig.[TJ The coordinates, sizes, and position angles 
for the apertures are also provided in Table [T] The orien- 
tation and the centre position of the aperture were chosen 
to derive an estimate of profiles vertical to the main axis 
with the mean close to the centre. 

For simplicity for all three filters we have chosen a pixel 
grid in the aperture with pixel size corresponding to that 
of the image at 250 /im (6"). The surface brightness at the 
pixel coordinates of the aperture was derived using the IDL 
cubic interpolation routine (parameter set to -0.5). As the 
image pixel sizes at 350 and 500 ^m are larger by default 
(10" and 14") the produced aperture maps at longer wave- 
lengths are oversampled in respect to the original images. 
The surface brightness profiles and uncertainties were de- 
rived as the mean value along the aperture width and its 
standard deviation. 



2.1. Filament a 

Filament a lies in an area located between the Cocoon 
Nebula to the west and the Northern Streamer to the east. 
It has a bone-like structure with a narrow filamentary struc- 
ture in the centre and wider extensions at both ends. The 
whole filament appears rather isolated with a well defined 
central part and an almost perfectly flat background. We 
have chosen a part with a symmetric surface brightness 
profile. 



2.2. Filament b 

Filament b is located north of the lower main part of 
the Northern Streamer which extends from east to west. 
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RA: 21h45m32.1s DEC: 47"38'23.9" 
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Fig. 1. SPIRE images taken with SPIRE at reference wavelength 250 (im of four individual regions in IC 5146 showing 
the filaments and the apertures used to analyze the surface brightness profiles. 



Table 1. Scan parameters 



61. 


GLON 


GLAT 


RA 


DEC 


e [°] 


length ['] 


width 


a 


93°51'10.79" 


-4°50' 22.48" 


21h48m32.7s 


47° 26' 23.0" 


+17.23 


3.0 


0.8 


b 


93°35'28.02" 


-4°21'38.94" 


21h45m32.1s 


47°38'17.9" 


-81.85 


2.4 


1.0 


c 


93°37'32.73" 


-4° 13' 29. 67" 


21h45ml0.1s 


47°45'51.3" 


-70.00 


3.0 


0.8 


d 


94° 27' 4.40" 


-5°39'13.64" 


21h54ml8.8s 


47°11'0.8" 


+40.51 


2.9 


0.8 



Interestingly, filament b lies almost perfectly orthogonal to 
the main structure. 



2.3. Filament c 

This filament app e ars a s number 5 in the study of 
I Arzoumanian et "ail (|201lD . It is located north-west from 
filament b and has a similar orientation. The filament was 
found to have a large column density and a rather small 
FWHM of ~ 0.06 pc (for a distant of 460 pc), in fact the 
smallest in their whole sample. 



The aperture contains additional emission in the east 
unrelated to the filament or the otherwise smooth back- 
ground. We modeled and subtracted this emission as de- 
scribed in App. [B] We restricted the fit by assuming a dust 
temperature of 15 K for this emission; we further fixed its 
central position in the direction along the width to be at 
the boundary of the aperture. 
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2.4. Filament d 

Filament d is a bright filament which extends to the east 
of the Cocoon Nebula. The bright emission and relatively 
warm temperature (Sect. 0]) indicates strong illumination 
by newly formed stars. This is also indicated by the warm 
dust temperature of 22.2 K derived for a filamentary struc- 
ture locate d close by to the south- west (filament 18 in 
the work of lArzoumanian et al.l (|201lD V The radiation field 
heating the dust grains in the filament might therefore be 
rather different in spectral shape and strength. It is likely 
also much less isotropic than in the case of the other three 
filaments. The incomplete sampling of the SED is more 
challenging because the temperature is derived on a less 
strongly curved SED. 

3. Model 

The surface brightness profile is in general a result of several 
different effects which are related to the intrinsic density 
profile, dust properties, illumination of the filament, and 
in general also temperature variations inside the structure. 
The method applied in this paper should be considered as a 
first approach, rather than a thorough analysis of the entire 
problem. In the following we describe the model assump- 
tions. 

3.1. Dust emissivity 

As mentioned in Sect. [H the surface brightness profiles 
in the SPIRE submm bands are less strongly affected by 
the decreasing temperature towards the centre of opaque 
clouds than are surface brightness profiles at the shorter 
wavelength PACS bands. The SPIRE observations are quite 
suitable for tracing the column density in the context of a 
single-temperature model, allowing us to ignore in this pilot 
study the radial temperature variation inside the cloud. 

With this common assumption that the emission from 
dust grains can be described by a single dust temperature, 
the emissivity per gram of matter (dust and gas) is given 

by 

Vu = B u (T dust )SK c ™ (j^) , (1) 

where B„(Td ust ) is the Planck function, Td us t the dust 
temperature, Kg™ the effective dust emission coefficient at 
Ao = c/vq = 250 /im where c is the velocity of light, and 
/3 em the corresponding effective exponent of the frequency 
dependence. 6 is the dust-to-gas ratio in mass for which 
we assumed 5 = 0.00588 as adopte d from a previo us study 
about the SED of condensed cores ()Fischerall2011|) . 

As the scattering at submm wavelength for interstellar 
dust grains is negligible the mean extinction is close to the 
mean absorption (or emission) properties. For single dust 
temperatures we would have k^" ~ where the extinc- 
tion coefficient is given by K^ xt = Kg xt [v / Vo)^"* . In the 
literature the product 6nf u in Eq. [T]is therefore replaced 
by the opacity SK CVLt . 

Here, we explicitly distinguish between the intrinsic 
dust properties which might even vary inside the cloud and 
the effective ones which correspond to an approximation 
where the dust emission peak is modeled by a simple mod- 
ified black-body function. 



Temperature variations inside the cloud will for example 
broaden the SED leading to an apparent lower /3 om . As 
the dust emission is dominated by the warm dust at the 
clouds outskirts the single dust temperature assigned to 
describe the spectral shape is too high and will lead to an 
underestimate of the dust mas s unless the effective dust 
emission coefficient is reduced (|Fischerall20TTl 12012ft . The 
effect can be a factor 2 to 3 depending on the central column 
density and the external pressure. 

In radiative t ransfe r ca lculations perfo rmed by 
iFischera k. Dopital (|2008t) and IFischeral (|2011[ ) the dust 
properties were chosen to be close to the properties in the 
diffuse phase with n c xt = 5.232 cm 2 g" 1 and /3 cxt = 2.024. 
For the effective emission behavior here we adopted a 
constant Kg m = 10 cm 2 g _1 and a spectral index /3 cm = 1.8. 
The higher constant was chosen to compensate for the 
effect of possible grain growth inside the cloud and the 
lower exponent for possible temperature variations. We 
will see that the adopted product Skq™ — 0.0588 cm 2 g _1 
is in broad agreement with the range of distances derived 
for the cloud complex (Sect. 15.1. ip . 

3.2. Mass surface density 

The filaments are assumed to be isothermal, self- 
gravitating, and in pr essure eq uilibrium with the ambient 
medium. As shown in baper II the shape of the density or 
pressure profile of an infinitely long cylinder can be char- 
acterized through the mass ratio / cy i = (M//)/(M//) max 
given as the ratio of the mass line density Mjl and the 
maximum possible mass line density given by 

IK 

= ~qi ( 2 ) 

where G is the gravitational constant. The constant K 
is given by K = kT cy \/ '(/ttmn) where k, /i, toh are the 
Boltzmann constant, the mean molecular weight, and the 
hydrogen mass. The gas is assumed to be molecular giving 
(i « 2.36. The effective temperature T cy \ contains the ther- 
mal and the turbulent motion of the gas. In the case of a 
non turbulent filament K — c 2 where c s is the sound sp eed. 

For T cy \ we assume 10 K. As discussed in Ipaper II this 
appears to be appropriate generally for dense cold struc- 
tures based on molecular line observations. We also found 
this model of pressurized isothermal self-gravitating fila- 
ments is in agreement with observations of the FWHM 
and central column density ATjj(O) of low mass filaments. 
This also indicates that the underlying hypothesis, that 
these filaments are in hydrodynamical equilibrium, is self- 
consistent, which motivates using this model in the ex- 
tended analysis in this paper. It would be ideal to have 
diagnostic molecular line observations to examine the tem- 
perature and dynamical state, as indeed there should be in 
some Gould Belt Survey fields, but to our knowledge there 
are none available for the filaments studied here. But be- 
cause T cy \ is an explicit scaling parameter of the model, 
the systematic dependence of the results on this choice can 
be examined (Sect. [5]). For 10 K, the maximum mass line 
density is 16.4 M /pc. 

The filament is assumed to be in pressure equilibrium 
with the surrounding medium with pressure p cx t- The mass 
ratio is related to the overpressure by p c /Pcxt = 1/(1 — /cyi) 2 
where p c is the central pressure in the cylinder. The mass 
ratio determines furthermore the shape of the mass surface 
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density while the external pressure its amplitude. The mass 
surface density at normalized impact param eter x = r / r cy \ 
where r cy i is the cylinder radius is given by (|paper J ) 



Poxt 



V8 



1 - /cyl 



4ttG 1 - / cyl 1 - / cyl (l - X 2 ) 

|v / /cyi(l-/c y i)(l-.x 2 ) 



1 - /cyl 



1 - /cyl(l - X 2 



x arctan < 



' /cyl(l-^ 2 ) 
1 - / cyl (l - X 2 ) 



(3) 



The corresponding column density of H nucleons is given 
by 



N H (x) = E m (x)(^to h ) \ 



(4) 



where /2 ~ 1.4. 

In the model the sky coordinates are related to the im- 
pact parameters using a stretch factor s = D/r cy \ so that 



S X (l? - 1? ) = »"A(gr] 



(5) 



where $0 is the central position of the filament. 

As noted in Ipaper J the mass surface density profile at 
given impact parameter does not depend on K but for given 
mass ratio and external pressure K affects the size of the 
cylinder. The radius for example is given by 



r cy i 



V8K 




(0) 



so that r cy i cx T cy \ / -^/Poxt for given f cy \. The radius is sym- 
metrical around f cy \ = 0.5 where the cylinder has the maxi- 
mum physical extension. The FWHM shows a qualitatively 
similar behavior as a function of the mass ratio. However, 
the maximum value appears at smaller f cy \ and at large 
/ cy i the FWHM decreases as (1 — / cy i) and is small relative 
to the physical radius. 

3.3. The SED 

The intrinsic surface brightness of a filament seen at incli- 
nation angle i at angle $ is given by 



1 



■£m(i%w 



(7) 



In our sample the profiles are affected by the finite res- 
olution. To allow a proper estimate of the intrinsic den- 
sity structure and the size of the filament we convolved 
the theoretical profile with the corresponding point spread 
function (psf) of the broad band filters. We used for sim- 
plicity the Gaussian approximation where the FWHM of 
the psf of the three SPIRE filters at the three reference 
wavelengths A r = 250 /xm, 350 /xm, and 500 /im is given 
by 17.6", 23.9", and 35.2", respectively (Table 5.2, SPIRE 
Observers' Manual, vers. 2.4, June 7, 2011). 

For the assumed model the flux density per unit length 
is simply given by 



Fl= I M = 



SM/l 
D cosi 



Vu = 



Sr] v 



D cos i 



7/ 



cyl" 



2K 



(8) 



where D is the distance to the filament. 

An important quantity in the interpretation of the dust 
emission is obtained by integratin g over frequen cy which 
gives the luminosity to mass ratio (|Fischerall201l| ) 



L/M 



2nD2F' 



5 4 (Kpf(T dust )) <7 SB IdustJ 



(9) 



(M /I) /cosi 

where ctsb is the Stefan-Boltzmann constant and where 

( Kpl (T dust )) = « ^— j ^_^r dust , (io) 

is the Planck averaged effective emission coefficient where 
a = 4 + f3 cm and where h is the Planck constant. We use 
this ratio to validate our assumption for the effective dust 
emission coefficient (Sect. [57 



3.4. Modeling profiles 

The surface brightness profiles for the model filament are 
determined by the dust temperature Td us tj the external 
pressure p ex t, the overpressure p c /p C xt or the mass ratio 
/cyi, the stretch factor s, and a position $ . It is assumed 
that the emission at reference frequency v r (r = 1,2,3) 
from the filament is superimposed on a smooth background 
of the surface brightness I^ rd (^) modeled through a poly- 
nomial function of first (filament a,b, and d) or second order 
(filament c). 

The profiles were modeled using a non-linear x 2_n t min- 
imizing 



3 

EE 



(ii) 



where I Vr {&i) is the observed surface brightness at position 
•di and J^(i9$) the corresponding theoretical surface bright- 
ness from the filament (see App. ITCJ corrected for spectral 
shape across the finite bandpass. In the modeling of the 
SED we accurately calculated the corresponding theoreti- 
cal fluxes for the broad band filters using the appropriate 
filter transmission curves. The color corrections are larger 
at longer wavelengths where the broad-band fluxes are more 
than ~ 10% higher than the corrected values. 

The dust temperature itself is sufficiently enough de- 
termined through the three broad-band filters. Otherwise, 
there are basically two main parameters determining the 
fit to the profiles: the mass ratio / cy i for the shape of the 
profile and the pressure which determines for given / cy i the 
amplitude. From the modeling standpoint we found that it 
is advantageous to perform the modeling in the parameter 
plane spanned by the central mass surface density Sm(0) 
(from Eq. [3]) and the mass ratio / cy i. 

4. Results 

Overall, we found a good agreement of the theoretical and 
the observed surface brightness profiles with a reduced 
Xred < 1 (Fig. [H in all four cases. The derived parame- 
ters for the four filaments are listed in Tab. [5] and Tab. [31 
The tem peratures for the fil a ments are close to those de- 
rived by lArzoumanian et al.l (1201 If) . The y are larger than 
the adopted T cy i. As discussed in baper 1 efficient coupling 
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Table 2. Fitted parameters using the model of a self- gravitating cylinder 



fil- Xrcd Nfrcc Tdust D COS i log 1Q / cyl 

[K] [pc] [log 10 Kcm- 3 ] 

no distance constraint 



a 
1) 
c 
d 


0.63 
0.29 
0.84 
0.16 


79 
61 
76 
76 


14.8 ±0.2 

13.9 ±0.5 
13.0 ±0.3 
19.9 ±0.7 


493 ± 36 
3771™ 
488 ± 31 
706 ± 70 


4.35 ±0.13 

°- Ui -0.45 
o ^7+0.36 
-3.47 

3.63 ±0.24 


0.68 ±0.05 
n 4Q+0-19 

u -^ y -0.43 

Q4+ 06 
0.81 ±0.05 








distance constrained (D = 


493 ± 36 pc) 




b 
d 


0.30 
0.24 


62 
77 


14.1 ±0.4 
18.9 ±0.5 


488 ± 32 
532 ± 31 


4.71 ±0.11 
4.16 ±0.16 


0.62 ±0.07 
0.70 ± 0.07 



Table 3. Derived parameters based on Tab. [2] 



fil. 


(M/0/cosi 


L/M 


FWHM 


FWHM cos i 


iV H (0)/cosi 




[M /pc] 


[Lq/Mq] 


['] 


[0.1 pc] 


[10 21 cm- 2 ] 


no distance constraint 


a 


11.0 ±0.8 


0.36 ±0.02 


0.63 ±0.03 


0.90 ±0.04 


10.4 ±0.5 


b 


8 1+ 3 ' 
8-J--7.0 


0.27 ±0.04 


0.54 ±0.09 


u -°°-0.49 


12.6 ± 1.8 


c 


1 >s 4+ - 9 

10 - 4 -0.5 


0.18 ±0.03 


0.35 ±0.02 


0.49 ± 0.05 


22.7 ±2.4 


d 


13.3 ±0.8 


2.18 ± 0.31 


0.61 ±0.03 


1.24 ±0.11 


8.5 ±0.8 


distance constrained (D = 493 ± 36 pc) 


b 


10.1 ± 1.2 


0.29 ± 0.03 


0.48 ± 0.03 


0.68 ±0.05 


12.9 ± 1.7 


fl 


11.4 ± 1.1 


1.61 ±0.18 


0.67 ± 0.03 


1.04 ±0.06 


9.2 ±0.9 



Table 4. Flux v r F l Vr [10- 10 W/m 2 ] 



fil. 



250^m 



350/im 



500^rn 



no distance constraint 

a 1.688 ±0.023 0.790 ± 0.003 0.264 ± 0.002 

b 1.301 ± 0.043 0.643 ± 0.007 0.222 ± 0.007 

c 1.441 ±0.103 0.776 ±0.045 0.283 ± 0.015 

d 4.352 ±0.104 1.556 ± 0.019 0.434 ± 0.010 

distance constrained 

b 1.324 ±0.027 0.645 ± 0.006 0.221 ±0.007 

d 4.150 ±0.066 1.543 ±0.020 0.441 ± 0.010 



is not expected in these filaments, especially in the out- 
skirts which dominate in determining Tdust • Because T cy i is 
a scaling parameter of the model, the choice does not affect 
the quality of fit of the model, but it systematically affects 
the parameters (Sect. [5] 

The fluxes of the broad band filters given in Tab. [4] are 
derived using Eq. [5] with the fitted results for / cy i, Dcosi, 
and Tdust (Fig. [3]). The three photometric images appear to 
be sufficiently accurate for a first analysis of the physical pa- 
rameters of the filaments despite the incomplete sampling 
of the SED. The largest systematic uncertainty is certainly 
expected for filament d. But even in this case the SED 
seems to be well determined. 

For comparison we also assumed a Gaussian profile 
which produced similar flux estimates. The derived pa- 
rameters for the dust temperature, size, and central col- 
umn density were all similar to the ones we obtained by 
a proper model of the emission profile. Larger deviations 
were found for filament c where the Gaussian approxima- 
tion leads to a systematically lower flux and an apparent 
larger size (FWHM). With this Gaussian approximation, 
however, we gain no information on D, or p e xt- In addition, 
in this approximation the mass ratio / cy i can only be deter- 
mined as the ratio of the derived mass line density and the 
maximum mass line density 2K/G, the former uncertain 
by the inclination angle and the distance D and the latter 



uncertain by T cy i. In the model fit presented here / cy i is 
determined from the shape of the profile. 

The result of the detailed model for the external pres- 
sure and the mass ratio / cy i is visualized in Fig. [4] We note 
that the regions of highest probability follow lines of con- 
stant central column density. The lines of constant central 
column densities are given by 



N H (0) 



at low mass ratios and by 



JVh(0) 



1 



Pcxt 



V8 



pm H V 47rG 1 - / cy i 2 



(12) 



(13) 



at high mass ratios / cy i >« 0.7. 

In the figure we added lines of constant size (FWHM). 
For given external pressure and cloud temperature the size 
varies as a function of mass ratio f cy \. At both low and 
high mass ratios the size approaches zero, at low mass ra- 
tios because of the vanishing mass and at large mass ratios 
because of compression. At / cy i 0.331 a cylinder has the 
maximum size at a given p cxt . Th e lines of constant FWHM 
are given through the equations (pa per J ) 



FWHM rz 
for / cy i <: 
FWHM R 



V3 



V8K 



0.1 and 



2V2273 



1- 



V8K 



V^GPc 



= (1 - /cyl), 



(14) 



(15) 



for /cyi >« 0.7. 

At high mass ratios for both constant FWHM and 
TVh (0) the functional dependence of p ex ± on / cy i is identical 
(p ex t oc (1 — / cy i) 2 ) and so the corresponding loci become 
parallel in Fig. |U 
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Fig. 2. Surface brightness profiles of the filaments at 250, 350, and 500 /zm, normalized by the central surface brightness 
of the intrinsic model, before convolution. The model is optimized so that after convolution it matches the data. Note that 
the multi wavelength data simultaneously fit the SED (Fig. [3]). The light grey shaded area shows the 1 sigma variations 
of the intrinsic surface brightness profiles of the model. The best fit model of the data has a peak less than unity because 
of convolution with the beam; the 1 sigma variations (most easily seen for filament b)) are shown by the dark areas. 
The data that were fit, corrected using the best fit for the background emission, are shown as points with error bars 
normalized to the best fit model. The adopted Gaussian approximations of the psf are given as dashed lines. North is to 
the right. Note that the scan lengths vary from filament to filament. 



The estimate of all three parameters (distance D, ex- 
ternal pressure p QX t, and mass ratio / cy i) requires accurate 
observations of the surface brightness profiles. To lower the 
uncertainty it is reasonable to use an additional constraint 
that the cloud distance should be roughly the same for all 
filaments. For filaments b and d we performed additional 
model fits where we constrained the distance using the es- 



timate derived for filament a; the distance is not fixed but 
allowed to vary within the la range. The corresponding fit- 
ted parameters are listed in Tabs. [2] and and the derived 
fluxes in Tab. |H As shown in Fig. [5] and Fig. [5] the ad- 
ditional constraint leads to tighter confidence intervals for 
the pressure and the mass ratio. 
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Fig. 3. Spectral energy distributions for the four filaments derived by a simultaneous fit of the profiles for unconstrained 
distances assuming for the dust emission spectra a modified black body. The theoretical flux densities of the SPIRE filters 
(grey filled circles) correspond to the best fit (black line). The vertical lines mark the uncertainties of the theoretical 
fluxes of the broad band filters and the horizontal lines the filter widths. The grey shaded areas show the la uncertainties 
of the SED. 



The filaments have deconvolved sizes (FWHM) on the 
order of 0.1 pc. All column densities are relatively high 
ranging from 7V H (0) ~8x 10 21 cm~ 2 to ~ 23 x 10 21 cm~ 2 . 
Filaments with higher central column density 7Vh(0) tend 
to have both a smaller FWHM and steeper density pro- 
file, characterized though a higher / cy i. This behavior is 
consistent with what is expected from the model of self- 
gravitating isothermal pressurized cylinders with mass ra- 
tios / cy i > 0.331. The iV H (0)-FWHM relation is discussed 
in Sect. O 

Comparing filaments a, b, and c we find a decrease 
of dust temperature with column density as expected 
for opaque filaments heated by the ambient radiation 
field. Strong radiative transfer effects are also indicated 
through the L/M ratio with values considerably lower than 
1.15 Lq/Mq, the typical r atio for diffuse d ust heated by the 
interstellar radiation field (|Fischerall2011l) . The filaments a, 
b, and c have cold dust temperatures below 15 K in agree- 
ment with what is expected for dust grains in the molecular 
phase. The warmer dust temperature and the higher L/M 
ratio for filament d indicate strong heating from recently 
formed stars in the Cocoon Nebula. The radiation field is 
probably therefore highly non isotropic. 

The four filaments indicate a trend in which more 
opaque filaments are also more compact as expected for fil- 



aments with strong self-gravity where FWHM oc 1/A r n(0) 
(Pig-IB. 

In the following we will describe the results for the four 
filaments in more detail. 



4.1. Filament a 

This filament provides the best estimate of the physical pa- 
rameters because of the well defined background and source 
emission as seen in Fig. [2J The surface brightness profile is 
well reproduced with the model. The only apparent devi- 
ation can be seen in the interior of the profile at 250 /jm 
which is asymmetric. This might be related to substructure 
below the resolution. Another possibility is the temperature 
variation caused by radiative effects and a non isotropic 
radiation field which is expected if the cloud complex as 
pointed out in Sect. [5] is located below the Galactic plane. 
In particular the UV and optical radiation could therefore 
be stronger at the surface facing the Galactic plane (right 
side in Fig. 0) leading to warmer dust temperatures. A sim- 
ilar asymmetry is seen in the surface brightness profiles of 
filament c. The origin of these asymmetries might be re- 
solved using PACS data. 
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Fig. 4. Regions of highest confidence in the f cy \-p e xt plane for the fit of the surface density profile at 250, 350, and 
500 fim with the model of an isothermal self-gravitating pressurized infinite cylinder (dark grey: 68%, medium grey: 90%, 
light grey: 99%). The cylinder is assumed to be seen edge-on. The corresponding overpressure p c /p ex t = (1 — /cyi)~ 2 is 
labeled on the right axis. The black curves give the column densities of neutral hydrogen through the centre, labelled 
with log 10 -/Vh(O) [10 21 cm -2 ]. The grey lines give the FWHM in parsec. 



4.2. Filament b 

The filament b shows the strongest variation of the surface 
brightness along the filament axis. A closer examination 
suggests that the variations are related to a condensation 
in the filament centre. Applying the method described in 
App. E] for the considered scan aperture we extracted the 
source emission. The source is elongated along the filament 
with axial ratio 2 and a length of 19.3" ± 0.7". The dust 
temperature with Tdust = 12.7 ± 0.3 is considerably colder 
than obtained for the filament. The mass is 1.01 ±0.07 M Q . 

Because of the larger uncertainties a precise determi- 
nation of all the parameters, distance, pressure, and mass 
ratio, is not possible. Fig.@]shows therefore a very extended 
confidence region. The distribution, however, is not random 
but follows a line of constant central column density. 

The fitted distance varies systematically with f cy \ and 
Pcxt- For given mass ratio f cy \ the distance decreases 
for higher pressure because the filaments are intrinsically 
smaller in higher pressure regions (Eq. [5]) . At low mass ra- 
tios (/ C yi < 0.5) the radius decreases for given pressure with 
/ cy i which leads to a systematic decrease of the distance 
with decreasing f cy \. Therefore, when we use the larger 
distance derived for filament a as a constraint we obtain 



more tightly constrained intrinsic profiles (Fig. [5]) and more 
tightly defined estimates of the gravitational state f cy \ and 
the external pressure (Fig. [5]), near the higher f cy \ and lower 
p C xt end of the original confidence range. 

4.3. Filament c 

Filament c is exceptional in several regards. It has the high- 
est central column density, the highest mass line density, 
the coldest dust temperature, and the smallest FWHM. 

The surface brightness profile indicates a filament with 
an extremely high mass ratio with / cy i > 0.8. The filament 
is therefore characterized by a high overpressure > 25 and 
by an intrinsic density profile at the outskirts steeper than 
r~ 2 . As we can see from Fig. [3J at large impact parameters 
the profile differs quite substantially from a simple Gaussian 
profile, showing wide tails. 

Even at the shortest wavelength the filament is not well 
resolved, with a size almost identical to the angular resolu- 
tion. Our model of the profile might therefore be affected 
to a larger degree by deviations of the true psf from the 
Gaussian function used here for sake of simplicity. 

An additional uncertainty is related to the source in the 
aperture that we subtracted using the method described in 
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Fig. 5. Same as Fig. [2] but for constrained distance 493 ± 36 pc. 




Fig. 6. Same as Fig. 2] but for constrained distance 493 ± 36 pc. 



App. [B] The excess seen at 1) « —0.75' in the background 
might be a residual of this source subtraction or a true 
structure. The interior profile shows a slight asymmetry 
with higher emission to the right. As in the case of fila- 
ment a that side faces towards the Galactic plane and the 
asymmetry might therefore be related to radiative transfer 
effects. 



5. Discussion 

As we have shown the dust emission from filaments in com- 
bination with a physical model can provide valuable infor- 
mation about the filament itself. The profile or the gravi- 
tational condition (overpressure) can be used to determine 
the pressure in the ambient medium and the distance to 
the filament. 



4.4. Filament d 



For this filament we obtain the largest distance. The value 
of Pext seems to be more in agreement with what is found 
for filament c. However, if we constrain the distance using 
the value derived for filament a we obtain for p cxt a value 
which is closer to that of filament a. The same applies for 
the mass ratio. 



5.1. Systematics 

The derived parameters for the distance and the pressure 
depend on the assumed inclination angle and the dust emis- 
sion coefficient. The distance further depends on the effec- 
tive temperature T cy \ which should ideally be determined 
in a separate observation using diagnostic molecular line 
observations tracing gas of mediu m and high densities a s 
used in studies of condensed cores (jRathbor ne et al.l [2008'). 
The model temperature T cy i could be higher than the ki- 
netic temperature if the filaments are also supported by 
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Fig. 7. Relation between the FWHM and the central col- 
umn density (Tab. [3l col. 4 and 6) for the four filaments 
in IC 5146 if the filaments are placed at the same distance 
D = 500 pc (black circles) closely consistent with the as- 
sumed effective emission coefficient Kg m = 10 cm 2 g _1 and 
an assumed dust-to-gas ratio S — 0.00588. For filaments b 
and d the values correspond to the fit with constrained dis- 
tance. The curves of the relation shown are der ived using 
the approximation (Eq. 32) provided in Ipaperl The grey 
shaded area shows the variation of the relation for a range 
of external pressures implied by the observed sizes and col- 
umn den sities of filaments in Polaris, IC 5146, and Aquila 
(jpaper J ). The external pressure is varied from 1.5 x 10 4 
to 5 x 10 4 K cm -3 where for given central column density 
the FWHM is larger for lower pressure. The given mass ra- 
tios on the upper axis correspond to an external pressure 
2 x 10 4 K cm -3 shown as the solid black curve. Also shown 
are the observed properties if the effective emission coeffi- 
cient is increased by a factor two (grey circles). The dashed 
curve shows the model relation forp e xt/fc = 5x 10 3 K cm -3 . 

magnet ic pressure (|Fieee fc Pudritdl2000allbh . However, in 
baper II we did not find any strong indication for an elevated 
T cy i that would support this possibility. We note, that in 
general, although neglected in this paper, the fit depends 
in addition on the exponent /3 Gm . 

A lower emission coefficient is, according to Eq. [T] and 
Eq. compensated by a higher mass surface density and 
therefore, because Sm oc -y/Pcxt (Eq. by a higher pres- 
sure. A similar effect is caused if the filament is seen at 
a projection angle i > (where i = refers to edge-on 
view). Thus, the derived pressure varies systematically as 

p cxt oc (cosi/^8 m ) 2 - 

This systematic dependence has a further implication 
for our distance estimate. According to Eq. [6] for a given 
mass ratio f cy \ the radius varies as r cy \ cx T cy \/ y?p cxt . Since 
D = s x r cy i our distances depend systematically as D oc 
T cy i5Kg m /cosi. The physical size (FWHM in pc) has the 
same dependence. 

The effect of doubling the emissivity on the relation 
between central column density and FWHM is shown in 
Fig. [7] The filaments have half the column density but are 
twice as extended. Wether the lower column density val- 
ues are consistent with the model might be verifiable using 
extinction measurements. 



For known D and T cy \ the model provides a method 
to study the dust properties and to estimate the effective 
emission coefficient in the dense molecular phase. The value 
is still uncertain because of the unknown inclination angle. 
On the other hand, the model would provide a method to 
estimate this angle if D, T cy i, Sk^" 1 , and /3 cm were known. 

5.1.1. Dust properties 

The dust properties, the dust composition and sizes, 
in the dense molecular phase are still rather uncertain. 
Theoretically, the opacity could, depending on composi- 
tion, be a factor ~ 10 h i gher than in the diffuse phase 
(jKrueeel fc Siebenmorgenl I1994D . According to theoreti- 
cal studies of grain coagulation the emission coefficient 
i n the FIR/submm regime should increase with density 
(jOssenkopf fc Henningj 119941) . Such a behavior is also in- 
dicated by recent studies of t he dust emission s pectrum 
for different column densities ()Martin et al.|[2012T ). In this 
case it might be expected that the dust emission coefficient 
for cold filaments would increase as a function of f cy \ if 
the bounding pressure is approximately the same. On the 
other hand, it needs to be considered that the effects of 
grain growth in dense clouds will be partly compensated 
by stronger temperature variations which lead to a reduced 
effective emission coefficient. The effect can be disentangled 
by using a more realistic model for the dust temperature 
inside the filament. 

Measurements of the opacity <5«; ext in the submm regime 
are generally based on correlations between the dust re- 
emission and the column density of the gas assuming = 
K| xt . The amount of gas and dust along the sight lines is 
derived using different tracers depending on the phase con- 
sidered. HI and CO observations are commonly used to 
trace the gas column density of the atomic and (up to a 
certain column density) molecular gas. Column densities in 
dense molecular gas can also be based on the reddening of 
background stars. In most recent publications the SED is 
modeled using a simple modified black body function as- 
suming a single temperature as in Eq. [TJ Therefore, the 
derived opacity needs to be considered as an approximate 
value. The correlation provides for given (3 cxt an estimate 
of the product Skq** or in units cm 2 /H — atom the opacity 
crff ' = I Amu. 

Based on recent studies of the optically thin medium 
< 7n xt lies in the range (0.6 — 1.6) x 10~ 25 cm 2 /H— atom 
()Planck Collaboration et ail 1201 If) . Our assumed value 
(a§ m = 1.38 x 10" 25 cm 2 /H - atom) lies at the high end of 
the derived range. In the dense molecular phase the effec- 
tive emission coefficient s eems to be a few time s higher than 
in the diffuse phase (see iMartin et al.l (|2012t ) for a recent 
summary and discussion). 

For known D and known T cy i our physical model of self- 
gravitating filaments provides an independent estimate of 
the dust opacity. As mentioned in Sect. [5] the previous dis- 
tance estimates vary within a factor of two and seem to 
depend on the applied method. The shortest distance of 
D = 460^go pc has been derived for the No rthern Streamer 
using star counts (jLada et al.l fl994. 1999), a value well in 
agreement with our results. Estimates based on the main 
sequence defined by B-type star s indicate for the cloud 
complex a distance around 1 kpc (jHarvev et al.ll2008h . The 
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model fit for filament a together with the range of distances 
suggests for the product Sk™ a physical range 

~ 0.055 cmV 1 < 5K c Q m ( — <~ 0.12 cx^g^tW) 
\ 10 Ky cos i 

where the lower value is related to the shortest distance 
460 pc. As T cy i in dense gas pro bably has a physical 
lower limit around 10 K ([paper If ) the model indicates 
an upper limit 5k™ 1 ~ 0.12 cm 2 g -1 or er™ 1 ~ 2.8 x 
10~ 25 cm 2 /H- atom. 

5.1.2. External pressure 

As we have shown in Ipaper J the physical size (FWHM) 
and the central column density iVn(0) of filaments of low to 
intermediate c olumn density in IC 5146, P olaris, and Aquila 
as derived by lArzoumanian et all ([20111 ) indicate T cy \ ~ 
10 K and a pressure range ~ 1.5 x 10 4 to 5 x 10 4 K cm -3 . 
For the same T cy \ and distance, we obtain similar values of 
p ex t here. However, it should be remarked that the column 
densities in the aforementioned work were based on Jk™ 1 = 
0.14 cm 2 g _1 whereas in the modeling here we used (5 k™ 1 = 
0.0588 cm 2 g- 1 ^ 

In a comp arable study about t he density structure of 
Bok-Globules iKandori et all ((2005) derived a range of ex- 
ternal pressure values ranging from 2.1 x 10 4 K cm -3 to 
1.8 x 10 5 K cm -3 . They agree within a factor of two with 
our study. 

Inspection of Figs. S] and E] reveals the possibility of a 
common p ex t/k ~ 2 x 10 4 K cm -3 , with filaments 'b' and 
'c' being on the high and low side, respectively. However, 
we feel that in the real ISM surrounding these different fil- 
aments the pressure might actually vary: some filaments 
might be embedded and therefore be surrounded by a 
medium with higher pressure. The model applied here and 
in other fields can then be used to understand how uniform 
the pressure is in the ISM within a cloud complex. 

A factor two higher emission coefficient consistent with 
the larger distance would imply a considerably lower exter- 
nal pressure for IC 5146 (p ext /k ~ 0.5x 10 4 K cm -3 ). As the 
larger distance also implies a large distance (—61 pc) below 
the Galactic mid-plane (Sect. [5]) a lower pressure would be 
also expected considering the pressure profile of the ISM 
(|Ferrierell200Tl ). 

5.2. Luminosity to mass ratio 

Radiative transfer calculations provide a further method to 
verify the assumption made for the effective dust emission 
behavior (k™ 1 and the exponent f3 cm ). Here, we used the 
L/M ratio to validate the choice of k™ 1 . 

Re levant studies can be found in the work of iFischeral 
(|2012f) . The estimates of the L/M ratio are based on 
radiative transfer calculations through pressurized self- 
gravitating isothermal cylinders and spheres which are illu- 
minated by the interstellar radiation field assumed to be 
isotropic. The calculations are based on dust properties 
close to the ones found in the diffuse interstellar medium. 
The intrinsic mean properties are identical with the ones 
mentioned in Sect. 13.11 

2 In Ipaper j we corrected N-q ( 0) for the wrong fi assumed in 
the work of lArzoumanian et all (|201ll ). 



Considered are clouds with central extinction values in 
the range from Ay = 0.1 mag up to Ay = 64 mag for 
two different external pressures: p c *t/k = 2 x 10 4 K/cm 3 
and p xt/fc = 10 6 K/cm 3 . The results for the high pressure 
regime do not seem to be of any relevance for the filaments 
analyzed in this work. 

For cylinders with a central extinction Ay = 0, 0.1, 
1, 2, 4, 8, and 16 mag the corresponding L/M ratios are 
1.15, 1.051, 0.835, 0.685, 0.530, 0.396, and 0.200 L Q /M Q , 
respectively. As expected, the ratio decreases towards fila- 
ments with higher column density. The radiation absorbed 
scales with the adopted interstellar radiation field. 

For the assumed effective emission coefficient the de- 
rived column densities for filaments a, b, and c cover a range 
from Ay ~ 5.6 mag to Ay ~ 12.1 mag assuming, as in the 
radiative transf er calculations, a total-to-selective extinc- 
tion Ry = 3.1 (|Fitzpatricklll999f) and a gas-to-dust rati o 
Nn/E(B-V) = 5.8x 10 21 cm 2 mag" 1 dBohlin et al.lll978h . 
The corresponding L/M values range from 0.363 ± 0.023 
to 0.180 ± 0.025 Lq/Mq. They are very close to the the- 
oretical values of the radiative transfer calculations. This 
is reassuring, but might well be fortuitous. The radiation 
field illuminating the filament as well as the absorption and 
scattering properties of the dust grains at UV and optical 
wavelengths might be different. In addition the extinction 
is uncertain by the unknown inclination angle. 

If we assume the larger cloud distance and therefore the 
higher effective emission coefficient the column density is 
lower by a factor two. The L/M ratio is (according to Eq.|H|) 
increased by the same factor. For filament a for example 
L/M = 0.72 ± 0.04 Lq/M g at Ay ~ 2.8 mag. In this case 
the observed ratio is ~ 50% higher than the corresponding 
theoretical value. However, this might be resolved given the 
uncertainties mentioned above. 

6. Summary 

We have analyzed the structure of four compact filaments in 
the cloud complex IC 5146 using recent Herschel observa- 
tions taken with the SPIRE instrument which provides the 
surface brightness at the reference wavelengths 250, 350, 
and 500 /im. Scans of the surface brightness profiles were 
derived based on four aperture windows positioned over 
parts of the filaments with well defined structure and back- 
ground. 

For predicting the filament emission we applied the 
model of isot hermal self- gravitating pressurized cylinders as 
described in Ipaper H . In this model the shape of the density 
profile is determined by the ratio / cy i = (M/l)/(M/l) max of 
the actual and the maximum possible mass line density. For 
simplicity we assumed the dust temperature to be constant 
throughout the filament. The dust emissivity is approxi- 
mated by a simple modified black-body spectrum with an 
effective dust emission coefficient approximated through a 
power law k™ 1 = k™ 1 ^ / Vo/ with Kq 111 = 10 cm 2 g _1 at 
250 /xm and f3 cm = 1.8. The dust-to-gas ratio in mass is 
assumed to be S = 0.00588. (The model depends on the 
product (5k™ 1 ) . The individual surface brightness profiles 
were modeled by taking the appropriate Gaussian approx- 
imations of the psf into account. 

Overall, we found good agreement between the 
modeled and the observed surface brightness profiles. 
In accordance with the physical model filaments with 
higher central column density Nn(Q) have a smaller 
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size (FWHM) and show a steeper density profile. 
Based on the best example we derived a cloud distance 
Dcosi = (493 ± 36)(fcg m /0.0588 erring" 1 ) (T cy i/10 K) pc 



(PextA)/ cos i 



and an external pressure 

10 4 - 35±013 (<5K§ m /0.0588 enr^g- 1 )- 2 K cm- 3 , where i 
is the projection angle of the filament with i = 
for an edge-on view. The previously estimated dis- 
tances of the cloud complex limit (5k™ 1 to the range 
- 0.055 cm^- 1 < <5Kg m (T cyl /10K)/cosi <~ 0.12 cnrV 1 - 
An upper limit is therefore Skq" 1 ^0.12 cm 2 g _1 . 

Considering the simplifications the results are encour- 
aging for further investigations. They should be based on 
a better resolution of the filaments (more nearby clouds), 
a complete coverage of the SED (including PACS), and a 
more realistic model of the temperature variation in the 
filament. This will provide a better estimate of the mass 
ratio / cy i, the external pressure and the distance, and pos- 
sibly allow investigation of the intrinsic dust properties. It 
might also probe to what extent the dust properties can be 
considered constant within a single dense filament. 
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Appendix A: Predicting surface brightnesses 

The filament is located in a distance where the resolution affects the 
observed profiles. To take this effect into account we used the Gaussian 
approximation of the psf of the SPIRE instrument. The theoretical 
surface brightness profile from the cylinder measured at a pixel reso- 
lution Ai? is then given by 



Appendix B: Source subtraction from the scan 

The structure of the ISM makes it difficult to find filaments on a flat 
or well defined background. We want to consider an aperture window 
containing, as in case of the aperture window chosen for filament c, not 
only the emission of the filament and a smoothly varying background 
emission but also of a marginally extended source. 

We modeled the source as an elliptical Gaussian function con- 
volved with the psf also assumed to be Gaussian. The major or mi- 
nor axis is further assumed to be aligned with the scan direction. The 
surface brightness from the source is therefore assumed to be 
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where F v is the flux density given by 
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0.0588 cm 2 g _1 and f3 = 1.8 and where uq = c/Ao 
where Aq = 250 fim. The standard deviations of the extra emission 
along the window width and the window length (x and y) are given 

by 
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where as and e<rg are the standard deviations of the source vertical 
and parallel to the scan direction and where tr ps f is the standard 
deviation of the psf . The source emission is described by the dust 
temperature T<j ust , the cloud mass M c \, the source size <rg, and the 
source aspect ratio e. 

To avoid any additional systematic uncertainties we made no fur- 
ther assumption for the surface brightness contributions from the 
background and the filament other than assuming that the surface 
brightness is constant along the scan width and only varies along the 
scan length. The surface brightness along the scan length is given by 
the mean 
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where l Vr is the theoretical monochromatic flux of the broad-band 
filter at reference frequency v r . In our calculations we made the con- 
volution at 10 time higher resolution than the pixel size of the scan. 

As the psf varies along the filter with wavelength and depends 
therefore on the dust emission spectrum it would be more accurate to 
derive the monochromatic flux for each filter after the convolutions. 
Here, we consider the effect to be small. A larger uncertainty might 
be related to use the Gaussian approximation instead of the actual 
psf. We note that our approach here is based on the assumption that 
the dust emission spectrum does not vary as a function of the impact 
parameter. This is certainly not true in the case of varying dust tem- 
perature that could exist because of radiative transfer effects which 
produce a declining dust temperature from the edge to the cloud cen- 
tre or because of a non isotropic radiation field. To provide the most 
accurate estimates of the surface brightness profiles one should derive 
the actual convolved profiles 7„ for a number of frequencies v 

which allow a determination of the monochromatic broad band fluxes. 

The parameters of the filaments were derived using a non-linear 
X 2 -fit where the partial derivatives for the parameters a-j were esti- 
mated in accordance with 
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quency v r and where N x is the number of pixel along the aperture 
width. 

The source parameters were derived from the residual 



^■Iv r {xi,yi) = I Vr (xi,yj) - (I Vr ) (Vj) 
by minimizing itcratively 
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